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MOTION OF AN AUTOGIRO ROTOR 
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SUMMARY 

An analysis is made of the influence on autogiro rotor 
characteristics of a periodic blade twist that varies with the 
azimuth position of the rotor blade and the results are com- 
pared with experimental data. The analysis expresses 
the influence of this type of twist upon the thrust, torque, 
and flapping motion of the rotor. The check against ex- 
perimental data shows that the periodic twist has a pro- 
nounced influence on the flapping motion and that this 
influence is accurately predicted by the analysis. The 
influence of the twist upon the thrust and torque could be 
demonstrated only indirectly, but its importance is indi- 
cated. 

INTRODUCTION 

The resultant of the air forces and the mass reactions 
on an autogiro rotor blade produces a couple tending 
to twist the blade unless the chordwise center of gravity 
of the blade and the center of pressure of the air forces 
are coincident. This fact has been known for some 
time, but it has not been generally realized that, except 
in particular cases, the resultant twist is periodic and is 
a function of the angular position of the blade in azi- 
muth. The periodic twist may be of a magnitude com- 
parable with or even exceeding the pitch setting of the 
rotor, which demonstrates the necessity of including it 
as a factor in the analysis of autogiro-rotor characteris- 
tics. 

It is the purpose of this paper to present an analysis 
of the periodic twist and to support the validity of the 
analysis by a comparison of predicted results with ex- 
perimental information obtained from flight tests of a 
direct-control wingless autogiro. The scope of the 
paper will be limited to a study of the influence of a 
known periodic twist upon the thrust, flapping motion, 
and torque of ° rotor; the effect of such a twist upon 
rotor vibrations and stability and the problem of pre- 
dicting the twist will be treated in a subsequent report. 

ANALYSIS 

It has been experimentally shown in previously 
unpublished data that, except in special cases, the air 



forces acting on a rotor blade cause a twist of the blade 
and a consequent change in the blade pitch angle. 
The twist is not constant but is a function of blade 
azimuth angle because of the variation of the air forces 
on the blade with this angle. In the following analysis 
the basic equations expressing the rotor characteristics 
will be generalized to include the factor of periodic 
twist. The notation of reference 1 will be used through- 
out this paper except for minor changes; for con- 
venience, the list of symbols is appended at the end of 
this section. 
The following additional notation is used: 

r=xR (1) 



dr=Rdx 



(2) 



The problem is now the solution of the equations for the 
autogiro rotor when the pitch angle 8 has the form 
0=0o+2#i +2*Q+2*i cos tf'+zih sin yff+xt % cos 2\p 

+X7b sin 2^+ . , . (3) 
where and 77* are coefficients descriptive of 0. 

It will be noted that an assumption is here made 
concerning the distribution of the twist along the 
radius; equation (3) shows that a linear distribution, 
starting from zero at the blade hub, has been used. 
Actually the twist reaches zero just outboard of the 
vertical pin and is not, in general, absolutely linear 
from there to the tip; the assumption used, however, 
does not introduce a serious error and is considered 
justified for its simplicity. 

It is obvious that the expressions of reference 1 for 
interference flow v, angle of attack a, blade flapping 
angle 0, and the dynamic equation of flapping are un- 
altered; they are 



2(\ 2 +m 2 )> 

0=<x o — a, cos ^—61 sin yf/—a 2 cos 2^ 
—b 2 sin 2^— .... 



(4) 
(5) 

(6) 
(7) 
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Using the notation of equation (1), the velocity com- 
ponents nt the rotor blade are: 



IIl=u T =x+n sin i> 



(8) 



^=« P =X+^i+(-M«o+^i+^) cos * 
+(-xa l +±iA^ sin *+(^i+2*& 2 ) cos 2* 
+Q M &! -2*0^ sin 2^+|m^ cos 3^+^*2 sin 3* 



(9) 



Also 



Ut »=^4-| m 2 +2mx sin *-km 2 cos 2^ 



(10) 



+(MX-a{x , -5M 2 ]-|^) ain* 

+ Oiara, + 2z 2 & 2 ) cos 2 ^ 2 Oo + »xb x - 22?ch) sin 2 $ 

-|^ 2 cos4^+|-A 2 sin4^ (11) 

The rotor thrust is calculated upon the assumptions 
that the elemental force on the blade lies in a pi** 16 
perpendicular to the blade-span axis and that the force 
depends only on the velocity in that plane; it is further 
assumed that u P is small compared with u T) so that the 
angle <p between u T and the resultant velocity (u T 2 + 
Up 2 )* may be equated to its sine and tangent. As an 
approximate allowance for tip loss, it will be assumed 
that the thrust becomes zero at a radius x= B where it 

is arbitrarily assumed that '5=1 — ^g- Then 

T^j^d^^lpcQWu/C.dx (12) 

It is further assumed that C L is a linear function of the 
blade-element angle of attack a n which is, of course, 
accurate below the stall; then 



C L —aa r 



(13) 
(14) 



Errors are introduced by the assumption that 
C L =aa T ; however, a graphical evaluation of the thrust 
made without this assumption and using a curve of C L 
against a derived from wind-tunnel tests disclosed that 



the error so introduced was negligible. Another error 
in the thrust expression exists where u T is negative; 
this error can be approximately nullified by the follow- 
ing correction. When u T is negative, the normal 
expression for the blade-element angle of attack must 
be altered to 

a/=-d-<p ' (15) 

and equation (15) must be used in the part of the disk 
bounded by x=-m sin * and and b ? * =ir and 
^=2t. The expression for the thrust is now, after 
substituting for C Li 0, and <p—u P lu Tj 

+xvi sin t+xti cos 2^+xij 2 sin 2*+ . . . )dx 

— ~pcan 2 R^ I ety I u T 2 (8*+x0i+x*o+ x ^ cos * 
+xm sin f+X€ 2 cos 2^+xi? 2 sin 2*+ ...)*! 

1 A /"*2» (*— m sin tf» 

-±pa&H?Z\ d*\ u T Urdx (16) 

r=i6 C p«ow(ix(B'+|/)+floQ^+^ 2 B-^M 3 ) 

+8^+^-^+^+^ 

It has already been shown (reference 1) that a H and 
b H are of the order m"; it can similarly be shown that e„ 
and rj n are of the order m r - 

The expression for thrust has been integrated upon 
the assumption, which experience has shown to be 
valid, that terms of higher order in m than the fourth 
are negligible. This same assumption will be used 
throughout the remainder of this analysis. 

The change in thrust caused by twist is 

(18) 

The thrust moment M T is, from reference 1, 
Mr = J[ B \ pcati 2 R* { 9u r * +u T u P }xdx, 

^2j[~" 3in ^pcafl 2 i?M^rH^p}^J r (19) 

where the second integral is added to thrust moment 
only in the interval *=t to *=2x. The integrated 
value of M r is 



Mr 
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12' 



- a,( J* - |m 2 5 2 ) - ^fi.B 3 + |^oB 4 
+ , I (^+J M 2 B 3 )-|^jsin * 
+ { - -J-mOoB 3 -0.03 5 M 'a 0 + 6,(^8' + ^B 2 ) 
-iMa 3 5 8 +«.(|B t +^S s )+|M'hB 4 ]co8 * 
+ [-]M 2 a»B 2 +^M 4 a 0 4-5M6 1 B 3 -ia 2 5 4 +iMt l B 4 
+n=(|B ! +|^)-|^jsin 2* 
+|-O.O53m 3 X-|M 2 0 o B 2 +^m 4 0o-!/W 

+ e2 (^+gM 2 B 3 )+j/«735 4 ]co8 2^ (20) 

Equation (7) can now be expanded into 

ISl\<h+3<h cos 2^+36 2 sin 2*) =M r -M (r (21) 

The coefficients of the flapping angle /J can now be 
obtained by substituting for M T in equation (21) and 
equating the coefficients of similar trigonometric func- 
tions; then, letting y=cpaR i lI l 

ao =I Y |i\F+0.080 M 3 X+i^4V5 , -gl' 4 ) 

+\™ Bi -h^\-m (22) 
«>=s^^i x ( B5 -H +e °(i 53+0106 ' l3 ) +<?lB4 

-\b i B s +^+ v j(l&+\sB^-\*B^ (23) 



(24) 



The expressions for a 2 and 6 2 must be expanded in 
powers of y. before solution is possible; it will be 
shown later that expansion to the order of n 2 is suf- 
ficient to express the thrust and torque to the order 
of m*- To this order the equations for a 2 and b 2 are 



6 ; 



16, 



b^B*=A-±a&+^&+\^+^M (25) 



1 m 



3 m ^4m' 



5 m 2 



4/i 5 m 



(26) 



Equations (25) and (20) may be solved for a 2 and 
b 2 after substituting for Oo, a u and bi to the appropri- 
ate order of n\ then 



a 2 = 



y 2 B*+ 144 



:K^ 7 -Sf) + <T +7 S?) 



+ ^12 + ^) + ^(l2 + ^) 

-i^+l^+TS'+IH (27) 



"7 2 B»-|- 14419 



15 



+ 5 7 m " + 30*. /J + 5 M 2 



_ 24 5.2 
"5yn 



(28) 



The changes in the expressions for the blade-motion 
coefficients arising from the twist are: 

^^yl^B'+l^yi^B'-^lp} (29) 

^-ff^-^+Ki'+i^-H (30) 
- 4 i(i B5+ r2^) + 2 1 4' lrf -l^) < 31 > 



iy 30 M 7 + 5 M 

(32) 



2 «t ; 



1 1l ; 



(33) 
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The air torque on the rotor blade is the sum of the 
accelerating torque arising from the lift elements and 
the decelerating torque arising from the drag elements. 
It has been assumed that the lift elements are zero 
between x—B and x=l; it is thought reasonable, 
however, to assume that the drag exists over the entire 
blade. An average value of the blade-element drag 
coefficient 5 will be used; this value is assumed con- 
stant with respect to the angle of attack. The torque 
Q, which must be zero by hypothesis, since the rotor is 
in a state of steady rotation, is then 

A flT TBI 

(2=0=^ I df J ~pcU 2 R 4 u T ^C L xdx 

-b- f*W Plpc&RWtodx (34) 
2ir Jo Jo * 

This expression after substitution for 0, C L , and p 
integrates into 

o= x {^-^)+x(^+| F MV.+^+i^ l ) 



+ 
+ 



^2(1^, 



(35) 



Equation (35) is a quadratic in X with the coefficients 
of the quadratic dependent upon y, B, 0i, and the 
« and t\ coefficients describing the periodic twist. The 
evaluation of equation (35) requires the substitution of 
known values of e and tj in the expressions for a n and 6„. 
The scope of this study includes only the prediction of 
the effect of a known twist upon rotor characteristics; 
consequently, the solution of equation (35) is possible 

when the drag term (l+p 2 —^ is known. 

Examination of equations (17) and (35) shows that 
the thrust and torque are expressed to the order m 4 if 
a 2 and 6 2 are expressed to the order /i 2 , inasmuch as 



and rf n are of the order m*. This last condition has been 
analytically proved but will not be included in this 
paper; the analysis of the e and jj coefficients of twist 
will be the subject of another paper. 

LIST OF SYMBOLS 

R, blade radius. 
6, number of blades. 
c, blade chord, feet. 
r } radius of blade element. 
S, r/R. 

$ 0 , blade pitch angle at hub, radians. 
0 t , difference between hub and tip pitch angles, 
radians. 

«„, coefficient of cos in expression for 0, radians. 
t;*, coefficient of sin in expression for 0, radians. 

0, instantaneous pitch angle, radians. 

5, mean profile-drag coefficient of rotor-blade airfoil 
section. 

\f/ t blade azimuth angle measured from down wind 

in direction of rotation, radians. 
v t rotor induced velocity. 

8, rotor angular velocity, d^/dt, radians per second. 
\UR, speed of axial flow through rotor, 
t&R, component of forward speed in plane of disk, 
equal to V cos a where V is forward speed, 
feet per second. 
f3 f blade flapping angle, radians. 
On, coefficient of cos n$ in expression for 0, radians. 
b n , coefficient of sin in expression for 0, radians. 
I i, mass moment of inertia of rotor blade about hor- 
izontal hinge, 
a, rotor angle of attack, radians. 
M T , thrust moment about horizontal hinge. 
M w , weight moment of blade about horizontal hinge. 
u T ttR, velocity component at blade element perpendic- 
ular to blade span and parallel to rotor disk. 
UpQR, velocity component at blade element perpendic- 
ular to blade span and to u T QR. 
T, rotor thrust. 
T 

Q, rotor torque. 

a, slope of curve of lift coefficient against angle of 
attack of blade airfoil section, in radian 
measure. 

a r , blade-element angle of attack, radians. 



cpaR* 



It 



f mass constant of rotor blade. 



c 

Z?=l— factor allowing for tip losses. 
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EXPERIMENT 

Flight tests were made of a Kellett KD-1 autogiro 
having the following characteristics: 

Gross weight, W 2,100 pounds. 

Rotor radius, R - 20.0 feet. 

Number of blades, b 3. 

Blade chord, c 1.00 foot. 

Blade weight, w b 61.5 pounds. 

Blade- weight moment, M w 482 pound-feet. 

Blade moment of inertia, 7 t _ . 175 slug-feet. 2 

be 

Rotor solidity, <r=^ft 0.0478. 

Blade mass constant, 

(sea level) y ==^-* 12.74. 

Blade airfoil section Gottingen 606. 

Pitch setting, $ Q (9 l = 0) 0.0960 radian. 

Airfoil section moment co- 
efficient, C ma c (about aero- 
dynamic center) - 0.056. 

Blade chordwise center-of- 
gravity location, c T (aft of 
aerodynamic center) 0.038 foot. 



240 



200 



&60 
a? 

Z/20 



80 



40 



































V 

0 








Rot(, 


y~ spe 


ed 




















































































Th 


rust c 


oeff/c 


tent 




















\ 





































.0048 



0044 



.o 

O040\ 
o 



5 



40 



60 80 
True air speed, m.p.h. 



100 



0036 



oas a 



Figure 1— Rotor speed and rotor thrust coefficient of KD-1 autogiro as measured 
In flight. 

The flight tests included the measurement of rotor speed 
as a function of air speed (fig. 1) from which, since the 
autogiro had no fixed wing, the thrust coefficient could 
be calculated. Simultaneous measurements of the 



blade flapping angle and twist were made with a motion- 
picture camera mounted on and turning with the rotor 
hub. The data for twist are shown in figure 2. The 
flapping-angle data are represented by the experimental 
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Figure 2.— Blade twist coefficients of KD-1 autogiro rotor as measured in Bight. 

points of figure 3. Both the flapping angle and the 
twist have been presented as the coefficients a n} b nf «„, 
and 77 „ of the expressions used in the previous section to 
represent 0 and 0. 

The effect of periodic twist upon the thrust coefficient 
C T was obtained by calculating the increment in C T 
caused by the periodic twist and deducting the incre- 
ment from the experimental value. The results are 
shown in figure 4. 

In order to check the derived expressions for the 
effect of periodic twist upon the flapping motion, the 
inflow factor X was calculated from the expression for 
the thrust (equation (17)) in which X was the only 
unknown; the calculation was made using the experi- 
mental values of the periodic twist, and it was also made 
on the assumption that all « and 17 coefficients except e 0 
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were zero. The two results are shown in figure 5. 
These values of X were then used in equations (22), (23), 
(24), (27), and (28) to calculate the blade flapping co- 
efficients both with and without the effect of the peri- 
odic twist; these results are shown in figure 3, together 
with the measured values of the blade-flapping-angle 
coefficients. Since the measured values of the periodic 
twist were obtained from targets placed at 3/4 R, they 
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Figure 3. -Blade flapping coefficients of KD-1 autogiro rotor as measured in flight 
and as calculated with and without periodic twist. 

were multiplied by 4/3 before insertion in the equations. 

The effect of the periodic twist in the torque equation 
was estimated in the following manner: Known values 
of X, a nf b Hf € nt rj n} B Qi and 9 t were substituted in equation 
(35) ; the resultant expression was used to evaluate the 

remaining unknown term ^(*+m 2 — £M 4 )- The ca *" 

culation was then repeated with the assumption that the 
periodic twist was zero and that X rather than the 5 term 



was unknown; the factors a nf b n , 9 0t 6 U and e 0 were 
assigned the same values in both calculations. The 
results of these calculations are showm in table I. 
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Figure 4.-Calcuhted periodic twist e'tcct on thrust cooilicient of KD-l autogiro. 
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Figure 5.— Inflow factor \ of KD-1 autogiro as calculated from thrust coefficient 
with and without periodic twist effect. 

DISCUSSION 

The influence of periodic blade twist on the rotor 
characteristics is illustrated in figures, 3, 4, and 5. 
The data in figure 3 afford convincing proof not only of 
the validitv of the twist analysis but also of the effect of 
periodic blade twist upon the rotor characteristics. 
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The data demonstrate that the type of twist developed 
in this rotor has a pronounced influence on the coning 
angle a 0 and on the flapping angle a x . The influence on 
the lag angle b x and on the second harmonics a? and b 2 
is considerably smaller. The agreement between the 
values of a Q and a x calculated from the expressions 
including the periodic twist and the experimental values 
is quite good. The calculated values of the lag angle 61 
are in radical disagreement with experiment; this same 
condition was encountered in previous work (reference 
1) and has been partially explained. In the reference 
it was shown that a variation of the rotor-induced 
velocity along the chord of the rotor disk had an appreci- 
able effect upon the variation of b x with p; an induced 
velocity increasing from the leading edge to the trailing 
edge increases b x . Since this type of asymmetry 
exists (reference 2) and varies inversely in magnitude 
with fi f the evaluation of its influence upon b x would 
improve the qualitative agreement between the calcu- 
lated and measured values. 

Figure 4 illustrates the magnitude of the periodic- 
twist contribution to the thrust coefficient. A different 
result was obtained in figure 5 by showing the difference 
in the values of X calculated when the periodic twist was 
considered and when it was neglected. 

In table I the calculated influence of periodic twist 
upon the torque equation and upon the resultant value 
of X is shown to be the least in magnitude of the effects 
studied. The effect is not, however, small enough to be 
neglected and would be an important factor if it were 
extrapolated to higher tip-speed ratios. 
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TABLE I. — EFFECT OF PERIODIC BLADE 
TWIST ON TORQUE EQUATION 



* 


X 

(experimental) 




X 

(calculated with- 
out periodic 
twist) 


0. 15 
.20 
.35 
.30 


a 0182 
.0186 
.0180 
.0180 


0.000700 
.000794 
.000845 
.000800 


0. 0177 
.0175 
.0171 
.0161 



1. The effect of periodic twist upon rotor-blade 
flapping coefficients is satisfactorily predicted by this 
analysis. 

2. The influence of periodic twist upon rotor char- 
acteristics as calculated from and checked with available 
data is an important factor in rotor analysis and can be 
adequately evaluated by the methods presented. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., January 28, 1937. 
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